Circle Packing for Origami Design Is Hard 

Erik D. Demaine* Sandor P. Fekete* Robert J. Lang* 
1 Introduction 

Over the last 20 years, the world of origami has been changed by the intro- 
duction of design algorithms that bear a close relationship to, if not outright 
ancestry from, computational geometry One of the first robust algorithms for 
origami design was the circle/river method (also called the tree method) de- 
veloped independently by Lang [THH] and Meguro [12l[l3]. This algorithm and 
its variants provide a systematic method for folding any structure that topo- 
logically resembles a graph theoretic weighted tree. Other algorithms followed, 
notably one by Tachi [15] that gives the crease pattern to fold an arbitrary 3D 
surface. 

Hopes of a general approach for efficiently solving all origami design prob- 
lems were dashed early on, when Bern and Hayes showed in 1996 that the 
general problem of crease assignment — given an arbitrary crease pattern, de- 
termine whether each fold is mountain or valley — was NP-complete [T]. In fact, 
they showed more: given a complete crease assignment, simply determining the 
stacking order of the layers of paper was also NP-complete. Fortunately, while 
crease assignment in the general case is hard, the crease patterns generated by 
the various design algorithms carry with them significant extra information as- 
sociated with each crease, enough extra information that the problem of crease 
assignment is typically only polynomial in difficulty. This is certainly the case 
for the tree method of design [3J. 

Designing a model using the tree method (or one of its variants) is a two-step 
process: the first step involves solving an optimization problem where one solves 
for certain key vertices of the crease pattern. The second step constructs creases 
following a geometric prescription and assigns their status as mountain, valley, 
or unfolded. The process of constructing the creases and assigning them is defi- 
nitely polynomial in complexity; but, up to now, the computational complexity 
of the optimization was not established. 

There were reasons for believing that the optimization was, in principle, 
computationally intractable. The conditions on the vertex coordinates in the 
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optimization can be expressed as a packing problem, in which the packing ob- 
jects are circles and "rivers, (which are curves of constant width) of varying 
size. It is known that many packing problems are, in fact, NP-hard, and our 
intuition suggested that this might be the case for the tree method optimization 
problem. 

In this paper, we show that this is, in fact, the case. The general tree method 
optimization problem is NP-hard. In the usual way with such problems, we show 
that any example of 3-Partition can be expressed as a tree method problem. 
At the same time, we show that deciding whether a given set of circles can be 
packed into a rectangle, an equilateral triangle, or a unit square are NP-hard 
problems, settling the complexity of these natural packing problems. On the 
positive side, we show that any set of circles of total area 1 can be packed into 
a square of edge length -j- = 2.2567 . . . 

2 Circle-River Design 

The basic circle- river method of origami has been described in we briefly 

recapitulate it here. As shown in Figure [TJ one is presented with a polygon P', 
which represents the paper to be folded, and an edge-weighted tree, T, which 
describes the topology of the desired folded shape. The design problem is to 
find the crease pattern that folds P' (or some convex subset) into an origami 
figure whose perpendicular projection has the topology of the desired tree T and 
whose edge lengths are proportional to the edge weights of T. The coefficient of 
proportionality m between the dimensions of the resulting folded form and the 
specified edge weights is called the scale of the crease pattern. The optimization 
form of the problem is to find the crease pattern that has the desired topology 
and that maximizes the scale m. 

Formally, the problem can be expressed as follows. There is a one-to-one 
correspondence between leaf nodes {rii} of the tree T and leaf vertices {vi} of 
the crease pattern whose projections map to the leaf nodes. We denote the 
edges of T by {e^} with edge weights w(e.j). For any two leaf nodes rii,nj G T, 
there is a unique path p$ j between them; this allows us to define the path length 
li j between them as 



We showed previously [3] that a necessary condition for the existence of a 
crease pattern with scale m was that for all leaf vertices, 



and subsequently, that with a few extra conditions, Equation [2] was sufficient 
for the existence of a full crease pattern (and we gave an algorithm for its 
construction). The largest possible crease pattern for a given polygon P' , then, 
can be found by solving the following problem: 




(1) 



\vi — Vj\ > ml 



(2) 



optimize m subject to 



{ 



\vi — vA > mli j for all i,j 
Vi G P' for all i 



(3) 
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Figure 1: Schematic of the problem, (a) P' is the paper to be folded, (b) T 
is an edge-weighted tree that describes the desired shape, (c) A solution to the 
optimization problem, showing creases and the ordering graph on the facets, 
(d) An x-ray view of the folded form, (e) A visual representation of the folded 
form. 
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Figure 2: A star tree and two possible solutions for the leaf vertices. Each 
solution corresponds to a packing of the circles centered on the leaf vertices. 

There is a simple physical picture of these conditions: if we surround each 
vertex by a circle whose radius is the scaled length of the edge incident to its 
corresponding leaf node and, for each branch edge of the tree, we insert into the 
crease pattern a curve of constant width (called a river) whose width is the scale 
length of the corresponding edge, then Equation [3] corresponds exactly to the 
problem of packing the circles and rivers in a non-overlapping way so that the 
centers of the circles are confined to the polygon P' and the incidences between 
touching circles and rivers match the incidences of their corresponding edges in 
the tree T. 

A special case arises when there are no rivers, i.e., in the case of a star tree 
with only a single branch node, as illustrated in Figure [2j In this case there are 
no rivers, and the optimization problem reduces to a single packing of circles, 
one for each leaf node, whose radius is given by the length of the edge incident 
to the corresponding node. 

Thus, several problems in origami design can be reduced to finding an op- 
timum packing of some number of circles of specified radii within a square (or 
other convex polygon). Several examples of such problems (and their solutions) 
are described in jlOj . 

We now show that this circle-packing problem is NP-complete. 

3 Packing and Complexity 

Problems of packing a given set of objects into a specific container appear in a 
large variety of applied and theoretical contexts. Many one-dimensional variants 
are known to be NP-complete, e.g., Bin Packing, where the objective is to pack 
a set of intervals of given lengths into as few unit-sized containers as possible. 
A special case of Bin Packing that is still NP-hard is 3-Partition, where an 
instance is given by 3n numbers xi with 1/4 < Xi < 1/2, and Y^i=i x i = n - 
Clearly, n unit-sized containers suffice for packing the object, iff there is a 
partition of the x% into n triples that each have combined weight 1; hence the 
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name 3-Partition. An important property of the problem is that it is strongly 
NP-complete: it remains hard even if there is only a constant number of different 
values Xi [5]. 

Like their one-dimensional counterparts, higher-dimensional packing prob- 
lems tend to be hard. Typically, the difficulty arises from complicated container 
shapes (e.g., a non-simple polygon to be filled with a large number of unit 
squares), or complicated objects (e.g., rectangles of many different sizes to be 
filled into a square, which is a generalization of Bin Packing.) This does not 
mean that packing simple objects into simple containers is necessarily easy: for 
some such problems it is not even known whether they belong to the class NP. 
One example is the problem Pallet Loading of deciding whether n rectangles 
of dimensions a x b can be packed into a larger rectangle of dimensions Ax B, 
for positive integers n, a, 6, A, B: it is open whether the existence of any feasible 
solution implies the existence of a packing that can be described in space poly- 
nomial in the input size log n + log a + log b + log A + log B, as the two different 
orientations of the small rectangles may give rise to complicated patterns. (See 
Problem #55 in The Open Problems Project, [3].) 

None of these difficulties arises when a limited number of simple shapes 
without rotation, in particular, different squares or circles are to be packed into 
a unit square. Leung et al. JT] managed to prove that the problem SQUARE 
Packing of deciding whether a given set of squares can be packed into a unit 
square is an NP-complete problem. Their proof is based on a reduction of the 
problem 3-Partition mentioned above: any 3-Partition instance LT3 P can 
be encoded as an instance LI sp of Square Packing, such that LI sp is solvable 
iff LT3 P is, and the encoding size of Yl sp is polynomial in the encoding size of 
LT3p. Membership in NP is not an issue, as coordinates of a feasible packing are 
integers of description size polynomial in the encoding size of LI sp . 

In the context of circle/river origami design, we are particularly interested 
in the problem of Circle Packing: given a set of n circles of a limited number 
of different sizes, decide whether they can be packed into a unit square. More 
precisely, we are interested in Circle Placement: given a set of n circles, 
place the circle centers on the paper, such that the overall circle layout is non- 
overlapping. Clearly, this feels closely related to Square Packing, so it is 
natural to suspect NP-completeness. However, when packing circles, another 
issue arises: tight packings may give rise to complicated coordinates. In fact, 
the minimum size C n of a C n x C n square that can accommodate n unit circles 
is only known for relatively moderate values of n; consequently, the membership 
of Circle Packing in NP is wide open. (At this point, n = 36 is the largest 
n for which the exact value of C n is known; see [14] for the current status of 
upper and lower bounds for n < 10,000.) 

Paradoxically, this additional difficulty has also constituted a major road- 
block for establishing NP-hardness of Circle Packing, which requires encod- 
ing desired combinatorial structures as appropriate packings: this is hard to do 
when little is known about the structure of optimal packings. 

The main result of this paper is to describe an NP-hardness proof of Circle 
Placement, based on a reduction of 3-Partition; it is straightforward to see 



5 



Figure 3: A symmetric 3-pocket with plug and shims. 



that this also implies NP-hardness of Circle Packing. In the following section, 
we will describe the key idea of using symmetric 3-pockets for this reduction: a 
triple of small "shim" circles C^, C'i 2 , Ci 3 and a medium-sized "plug" circles 
can be packed into such a pocket, iff the corresponding triple of numbers x^, 
Xi 2 , Xi 3 add up to at most 1. In the following sections, we show how symmetric 
3-pockets can be forced for triangular paper (Section [5]), for rectangular paper 
(Section O, and for square paper (Section [7]). The technical details for the 
proof of NP-hardness are wrapped up in Section [8] and Section [9j in which we 
sketch additional aspects of filling undesired holes in the resulting packings, 
approximating the involved irrational coordinates, and the polynomial size of 
the overall construction. On the positive side, we show in Section [TU] that circle 
packing becomes a lot easier if one is willing to compromise on the size of the 
piece of paper: we prove that any given set of circles of total area at most 1 can 
easily and recursively be packed into a square of edge length -4= = 2.2567 . . . 

4 Symmetric 3-Pockets 

Our reduction is based on the simple construction shown in Figure [3] It consists 
of a symmetric 3-pocket as the container, which is the area bounded by three 
congruent touching circles. 

Into each pocket, we pack a medium-sized circle (called a plug) that fits 
into the center, and three small identical circles (called shims) that fit into the 
three corners left by the plug. Straightforward trigonometry (or use of Proposi- 
tion !8.1[) shows that for a pocket formed by three unit circles, the corresponding 
size is 2/v3 — 1 = 0.1547... for the plug; the value for the shims works out to 
1/(5 + V3 + 2^7 + 4^3 = 0.07044.... 

Clearly, this packing is unique, and the basic layout of the solution does not 
change when the plug is reduced in size by a tiny amount, say, s = 1/N for a 
suitably big N. while each shim is increased by a corresponding amount that 
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keeps the overall packing tight. This results in a radius of r p for each plug, and 
a radius of r s for each shim. 

Now consider the numbers Xi for i = 1, . . . , 3n, constituting an instance of 
3-Partition. We get a feasible partition iff all triples (ii, 12,13) are feasible, 
i.e., Sj=i x ij = 1- By introducing x\ = 1/3 — Xi and using X2»=i x i = n > ^ is 
easy to see that a partition is feasible iff Y^j=i x \j — f° r au triples (ii, «2, *3)- 
Note that a 3-Partition instance involves only a constant number of different 
sizes, so there is some 5 > such that any infcasiblc triple (ii, 12, 43) incurs 
Sj=i x 'i — fi- By picking N large enough, we may assume 6 > e. 

As a next step, map each xi to a slightly modified shim Si by picking the 
shim radius to be r, = r s — x^/N 2 . We will make use of the following elementary 
lemma; see Figured] 

Lemma 4.1 Consider an equilateral triangle A = i>2> ^3) bounded by the 
lines £1,12,(3 through the triangle edges e\ 1 e2,ej,. For an arbitrary point p, let 
dj be the distance of p from ij. Define yj = dj, if p is on the same side of £j as 
A, and yj = —dj if p is separated from A by £j. Then Uj * s independent 

of the position of p. 

Proof: Consider the three triangles (v±,V2,p), (1)2,1)3, p), (i>3,vi,p). Their 
areas are d.3/2, di/2, g?2/2, hence j/i/2 + 2/2/2 + 2/3/2 is always equal to the area 
of A, i.e., a constant. □ 

The crucial argument for our reduction is the following. 
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Figure 5: Changing the size of a shim. 



Lemma 4.2 A set of three shims S^, Si 2 , Si 3 and a plug P of radius r p can 
be packed into a 3-pocket, iff X^j=i — 0; i-S-, if (il, i%, 23) is feasible. 

Proof: Refer to Figure [5] Let c be the center point of the pocket. For 
each of the three corners of the pocket, consider the two tangents and T 2 
between an unmodified shim of radius r s and the touching pocket boundary; 
let 2<f> S]0,7r[ be the angle enclosed by those two tangents. (If the pocket was 
an equilateral triangle, we would get <f> = 7r/6; the exact value for pockets with 
circular boundaries can be computed, but the exact value does not matter.) 

Now consider the shim motion arising by modifying r s by x\ . /N 2 , while 
keeping the shim tightly wedged into the corner. This moves its center point 
Pi, along the bisector b between T}. and T 2 . Let c = 1/ sin (f>. Considering 
the first-order expansion of the shim motion, we conclude that p,-. moves by 
c x x'ijN 2 + Q(l/N 4 ) along b, to a position qj. 

Finally, refer to Figure [S] and consider the possible placement of a plug after 
placing the modified shims , Si 2 , Si 3 into the the three corners; this requires 
finding a point within the pocket that has distance at least r p + r s — x\ . /N 2 
from each pi j . For this purpose, consider the circle Ci j of radius r p + r s — x\ . /N 2 
around each pi . . As shown in Figure [7J let U . be the tangent to Cj at the point 
closest to c; let di j be the distance of c to t^ . If we define yi . = di. for c is 
outside of Cj. , and = —d^ for c is inside of Cj. , then y^ = — ((c + 1) x 
x'iJN 2 + 9(1/7V 4 )). Now consider the set A of points separated by from 
Pi t , U 2 from pi 2 , ti 3 from pi 3 . Making use of Lemma I4.1[ we conclude that A is 



Figure 6: Finding a feasible placement for the plug. 




Figure 7: The existence of a feasible placement for the plug depends on the 
sum of distances of c from the sides of A. (Distances are not drawn to scale so 
that circles and tangents can be distinguished; in reality, they are much closer.) 



9 



Figure 8: The unique packing of (fc + 2)(fc + l)/2 unit circles into an equilateral 
triangle of edge length 2k leaves k 2 identical symmetric 3-pockets. 

a nonempty isosceles triangle, iff YlVij — 0> i- e -i ^ ^x^./N 2 — 9(l/iV 4 ) < 0. 
Given that J2 x ij > implies Yl x 'i — $ > V-^> we conclude that ^x^. < 
implies the existence of a feasible packing. 

Conversely, consider > 0- Given that each has a distance 0(1/7V 2 ) 

from c, we observe that the corners of the triangle formed by , ti 2 , U 3 are 
within 0(1/A f4 ) from Ci 11 Ci 2 , Ci 3 . However, because of Yl x i — $ > 1/-^) 
we conclude that any point of A is at least 0(1/7V 3 ) from being feasible. This 
implies that there is no feasible placement for the plug, concluding the proof. 

□ 

5 Triangular Paper 

For making use of Lemma 14.21 and completing the reduction, we need to define 
a set of circles (called rocks) that can only be packed in a way that results in a 
suitable number of 3-pockcts. In the case of triangular paper, this is relatively 
easy by making use of a result by Graham [5] . 

Proposition 5.1 An equilateral triangle of edge length 2k has a unique packing 
o/(fc + 2)(fc+l)/2 unit circles; this uses a hexagonal grid pattern, placing circles 
on the corners of the triangle. 

This creates 5^v =1 (2i — 1) = k 2 symmetric 3-pockets. After handling some 
issues of accuracy and approximation (which are discussed in Section [9]), we get 
the desired result. 

Theorem 5.2 Circle/river origami design for triangular paper is NP-hard. 
As a corollary, we get 
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Figure 9: Packing 2k unit circles into a rectangle of dimensions 2k - 1 and 
leaves 2k — 2 identical symmetric 3-pockets. 

Corollary 5.3 It is NP-hard to decide whether a given set of circles can be 
packed into an equilateral triangle. 

6 Rectangular Paper 

Similar to triangular paper, it is easy to force a suitable number of symmetric 3- 
pockets for the case of rectangular paper, see Figured] disregarding symmetries, 
2k unit circles can only be packed into an 2k — 1 by y/3 rectangle in the manner 
shown. This creates 2k — 2 symmetric 3-pockets, which can be used for the 
hardness proof. 

Because the input created for encoding an instance n 3p of 3-Partition 
needs to be a set of rationals whose size is bounded by a polynomial in the 
encoding size of n 3p , the irrational numbers needs to be suitably approximated 
without compromising the overall structure. This will be discussed in Section^ 
As a consequence, we get 

Theorem 6.1 Circle/river origami design for rectangular paper is NP-hard. 

This yields the following easy corollary. 

Corollary 6.2 It is NP-hard to decide whether a given set of circles can be 
packed into a given rectangle. 

7 Square Paper 

Setting up a sufficient number of symmetric 3-pockets for square paper is slightly 
trickier: there is no infinite family of positive integers n, for which the optimal 
patterns of packing n unit circles into a minimum-size square are known. As a 
consequence, we make use of a different construction; without loss of generality, 
our piece of paper is a unit square. 

As a first step, we use four large circles of radius 1/2, creating a symmetric 
4-pocket, as shown in Figure \W\ Now a circle of radius v/ ^~ 1 has a unique 
feasible placement in the center of the pocket, leaving four smaller auxiliary 
pockets, as shown. 
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Figure 10: A gadget for creating identical triangular pockets: the shown set 
of 13 circles has a unique packing into a symmetric 4-pocket. This creates four 
smaller symmetric 3-pockets, indicated by arrows. 

Now we use 12 identical "plug" circles and four slightly smaller "fixation" 
circles, such that three plugs and one shim have a tight packing as shown in 
the figure. For these it is not hard to argue that not more than three plugs 
fit into an auxiliary pocket, ensuring that precisely three must be placed into 
each pocket. Moreover, it can be shown that at most one additional shim 
can be packed along with the three plugs; this admits precisely the packing 
shown in the figure, creating a symmetric 3-pocket in each auxiliary pocket. In 
addition, we get a number of undcsired asymmetric pockets, which must be used 
for accommodating appropriate sets of "filling" circles, leaving only small gaps 
that cannot be used for packing the circles that are relevant for the reduction. 

As shown in FigureJTTJ we can use a similar auxiliary construction (consisting 
of 13 circles) for the 3-pockets in a recursive manner in order to replace each 
symmetric 3-pocket by three smaller symmetric 3-pockets. The argument is 
analogous to the one for 4-pockets. Again, additional filling circles are used; 
these do not compromise the overall structure of the packing, as the overall 
argument holds 

Theorem 7.1 Circle/river origami design for square paper is NP-hard. 

This yields the following easy corollary. 

Corollary 7.2 It is NP-hard to decide whether a given set of circles can be 
packed into a given square. 
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Figure 11: A gadget for creating multiple identical triangular pockets: the 
shown set of 13 circles has a unique packing into a symmetric 3-pockct. This 
creates three smaller symmetric 3-pockets, indicated by arrows. 

8 Filling Gaps 

The above constructions create a number of additional gaps in the form of 
asymmetric 3-pockets. Each is bounded by three touching circles, say, of radius 
r i: r 2, r 3 . By adding appropriate sets of "filler" circles that precisely fit into 
these pockets, we can ensure that they cannot be exploited for sidestepping the 
desired packing structure of the reduction. Computing the necessary radii can 
simply be done by using the following formula. 

Proposition 8.1 The radius r of a largest circle inscribed into a pocket formed 
by three mutually touching circles with radii r±, r 2; r 3 satisfies 

l/r = 1/n + l/r 2 + l/r 3 + 2 v /l/r 1 r 2 + l/nr 3 + l/r 2 r 3 

Note that the resulting r is smaller than the smallest n, and at least a 
factor of 3 smaller than the largest of the circles. Therefore, computing the 
filler circles by decreasing magnitude guarantees that all gaps are filled in the 
desired fashion, and that only a polynomial number of such circles is needed. 

9 Encoding the Input 

In order to complete our NP-hardness proof for Circle Placement, we still 
need to ensure that the description size of the resulting Circle Placement 
instance is polynomial in the size of the input for the original 3-Partition 
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Figure 12: A quad-tree packing guarantees that any set of circles of total area 
at most 1 can be packed into a square of edge length 7 = 4/y / 7r = 2.2567... 

instance. It is easy to see from the above that the total number of circles 
remains polynomial. This leaves the issue of encoding the radii themselves: 
if we insist on tightness of all packings, we get irrational numbers that can be 
described as nested square roots. As described in SectionHJ the key mechanisms 
of our construction still work if we use a sufficiently close approximation. This 
allows to use sufficiently tight approximations of the involved square roots in 
other parts of the construction, provided the involved computations are fast 
and easy to carry out. For our purposes, even Heron's quadratically converging 
method (which doubles the number of correct digits in each simple iteration 
step) suffices. 

10 A Positive Result 

Our NP-hardness results imply that there is little hope for a polynomial-time 
algorithms that computes the smallest possible triangle, rectangle or square for 
placing or packing a given set of circles. However, it is possible to guarantee 
the existence of a feasible solution, if one is willing to use larger paper. In fact, 
we show that a square of edge length 4/y / 7r = 2.2567... suffices for packing any 
set of circles that have total area 1 . 

Theorem 10.1 Consider a set S of circles of total area 1, and a square S of 
edge length A/^/tt. Then S can be packed into S . 
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Figure 13: A lower-bound example for packing circles: two circles of area 1/2 
require a square of edge length at least 1.362... 

Proof: Refer to Figure fT2l For each circle Ci of radius r^, let rii be chosen 
such that 2 "l+i < r i — th- Hence, replacing each Ci by a square Si of size 
increases the edge length by a factor of at most 7 = 4/ ^/tt. Now a recursive 
subdivision of S into sub-squares of progressively smaller size can be used to 
pack all squares Si, showing that all circles d can be packed. 

11 Conclusions 

In this paper, we have proven that even the aspect of circle packing in circle/river 
origami design is NP-hard. On the positive side, we showed that the size of a 
smallest sufficient square for accommodating a given set of circles can easily be 
approximated within a factor 2.2567... A number of interesting open questions 
remain: 

• Our 2.2567-approximation is quite simple. The performance guarantee is 
based on a simple area argument. This gives rise to the following question: 
what is the smallest square that suffices for packing any set of circles 
of total area 1? We believe the worst-case may very well be shown in 
Figure [T3J which yields a lower bound of (1 + y/2)/y/n = 1.362 . . . There 
are ways to improve the upper bound; at this point, we can establish 
2\/2/Vtt = 1.5957... [2]. 

• The same question can be posed for placing circles instead of packing 
them. 

• The approximation of circle packing does not produce a "clustered" layout 
as required by circle/river origami design, where objects that are close in 
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the hierarchy should be place in close vicinity. In the absence of rivers, 
we can reproduce the quad-tree packing in this context by making use of 
a space-filling curve. 

• On the other hand, we do not know yet how to approximate the necessary 
paper size in the presence of rivers of positive width. 
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